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DETERMINATION  OF  SOIL  PROPERTIES 


THROUGH  GROUND  MOTION  ANALYSIS 

John  Frye 
& 

Norman  Llpner 
S  June,  1981 
ABSTRACT 


A  method  of  calculating  In  situ  one  dimensional  stress-strain  soil 
properties  from  vertical  ground  motion  Is  presented.  The  method  relies 
on  the  fact  that  superselsmlc  air  blast  ground  surface  loadings  produce 
ground  motions  that  are  very  nearly  vertical  and  one  dimensional  In 
character.  Therefore  the  equations  of  motion  that  govern  the  response  are 
simple  and  may  be  Integrated  to  obtain  one  dimensional  stress-strain 
relations.  Thus,  results  from  tests  that  ’incorporate  superselsmlc  air 
blast  surface  loading  and  sensors  to  measure  vertical  motion  at  various 
depths  In  the  soil  can  be  used  to  calculate  soil  stress-strain  properties 
directly.  The  method  accounts  for  multiple  records  at  a  given  depth  and 
features  techniques  for  characterizing  response  histories  and  Interpolating 
velocities  at  depths  between  those  where  measurements  have  been  made.  As 
an  example,  for  the  DISC  HEST  Test  I  event,  conducted  In  Ralston  Valley, 
Nevada  as  part  of  the  MX  development  program,  the  site  properties  are 

computed  based  on  the  free  field  data. 

v 


DETERMINATION  OF  SOIL  PROPERTIES  THROUGH  GROUND  MOTION  ANALYSIS 
J.  W.  Frye1  and  N.  Llpner2!  M.  ASCE 

INTRODUCTION 

The  ground  response  to  overhead  highly  superseismlc  alrblast  loading 
(airblast  shock  speed  faster  than  ground  shock  speeds)  Is  nearly  one-dimensional 
uniaxial  strain  and  the  motions  are  nearly  vertical.  Soil  properties  for 
prediction  of  these  motions  have  typically  been  determined  from  dynamic 
uniaxial  strain  laboratory  tests.  However,  the  process  of  extracting  soil 
samples  from  the  field  can  disturb  the  material  and,  as  a  result,  the  laboratory 
properties  could  be  different  from  the  in  situ  material  behavior. 

An  approach  to  determine  the  in  situ  properties  is  to  perform  a  field 
test  where  the  ground  surface  Is  loaded  by  superseismlc  airblast.  Data  from 
sensors  that  measure  vertical  motion  at  various  depths  in  the  ground  could 
then  be  used  to  calculate  the  uniaxial  strain  properties  of  the  soil  by  use 
of  the  one-dimensional  equation  of  motion. 

One  type  of  surface  loading  that  has  been  used  to  obtain  a  one-dimensional 
response  is  the  DISC  (Dynamic  In  Situ  Compressibility)  HEST  (High  Explosive 
Simulation  Technique)  test  shown  in  Figure  1.  This  test  employs  a  circular 
region  of  explosives  that  is  center  detonated.  The  detonation  propagates 
outward  fast  enough  that  the  early-time  response  to  peak  velocity  is  essentially 
one-dimensional  within  some  region  under  the  loaded  area  that  Is  governed  by 
the  disc  radius  and  the  soil  properties.  Because  of  the  finite  propagation 
velocity,  time  at  any  range  from  the  centerline  is  measured  with  respect  to 
the  arrival  of  the  overhead  airblast  at  that  range. 

1.  Member  of  the  Technical  Staff,  Hardness  and  Survivability  Laboratory, 

TRM,  Redondo  Beach,  CA  90278 

2.  Department  Head,  Hardness  and  Survivability  Laboratory, 

TRM,  Redondo  Beach,  CA  90278 
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mWALYSIS  formulation 


The  one-dimensional  equation  of  motion  relates  the  vertical  normal 
stress  gradient  to  the  acceleration  of  the  soil. 
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Here  z  Is  the  vertical  coordinate  (Fig.  1),  o  is  the  normal  stress  in  the  s 
direction  (tensile  stress  is  positive),  and  p  Is  the  soil  density. 

Equation  (1)  can  be  integrated  with  respect  to  depth  from  the  ground 
surface  to  a  depth  z  to  obtain  the  following  equation  for  stress: 


-p(t)  p&z<fe 


(2) 


where  the  constant  of  integration,  p(t) ,  is  the  surface  pressure-time 
history,  a  boundary  condition  of  the  problem,  and  the  first  time  deriva¬ 
tive  of  velocity,  u  .  replaces  the  second  time  derivative  of  displacement. 
The  one-dimensional  strain,  e.  is  the  derivative  of  the  vertical 

Z 


displacement  with  respect  to  depth. 
3  «_ 


(3) 


Z  32 

Taking  the  time  derivative  of  Equation  (3)  provides  the  following  relation 


(4) 


for  the  strain  rate  e  : 
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and  integrating  Equation  (4)  gives  the  following  relation  for  strain  In  terms 
of  velocity: 


(S) 


The  constant  of  Integration  Is  zero  because  the  strain  in  the  soil  Is 


measured  with  respect  to  the  geostatic  strain  and ,  therefore.  Is  zero  at 
time  zero. 

From  Equations  (2)  and  (5)  It  is  seen  that  If  the  vertical  component 
of  velocity  is  defined  with  respect  to  time  and  depth,  then  stress  and 
strain  may  be  directly  calculated  from  derivatives  of  the  velocity  by 
performing  Integrations  with  respect  to  depth  and  time.  The  pressure 
time  history  at  the  surface  and  estimates  of  the  soil  density  are  also 
required  by  Equation  (2). 

Evaluation  of  the  Integrals  of  Equations  (2)  and  (5)  requires 
knowledge  of  the  velocity  field  for  the  complete  space-time  region  of 
interest.  In  a  test,  motion  sensors  at  only  a  limited  number  of  depths  can 
be  implanted  because  of  cost,  as  well  as  physical  constraints.  Therefore,  a 
method  must  be  developed  to  interpolate  from  the  data  available  for  a  limited 
number  of  depths  to  the  velocity  history  at  any  depth. 

In  many  test  events,  more  than  one  record  is  available  at  some  depths, 
so  that  if  one  sensor  Is  faulty,  all  of  the  velocity-time  information 
concerning  that  particular  depth  Is  not  lost.  Data  records  available  for 
a  particular  depth  will  vary  from  one  another  due  to  a  number  of  reasons, 
such  as  variation  from  one  location  to  another  of  soil  properties  and  surface 
pressure-time  histories.  In  examining  the  velocity-time  records  taken  at 
a  particular  depth,  it  is  not  always  obvious  that  one  particular  record  Is 
the  most  accurate  and  representative  of  all.  Thus,  some  method  of  Including 
all  acceptable  records  taken  at  a  particular  depth  must  be  employed  in 
defining  the  velocity-time  history  to  be  used  In  the  interpolation  process. 
Clearly  some  records  that  show  anomalous  results,  not  representative  of 
motions  that  are  physically  plausable,  should  be  excluded  from  consideration 
Such  records  might  be  taken  from  sensors  that  either  were  poorly  installed. 


were  destroyed  by  the  shock  loading,  or  produced  records  that  were 
extremely  noisy. 

The  method  chosen  for  interpolating  the  measured  velocity  records  uses 
a  transformed  coordinate  system  with  a  time-like  coordinate  a  that  remains 
constant  with  depth  along  specific  characteristics  of  the  velocity  time 
history.  At  the  beginning  of  motion,  a  is  always  taken  as  zero;  at  peak 
velocity,  a  equals  1.0;  and  at  the  end  of  the  velocity  record,  8  assumes 
some  large  value  such  as  10.  Other  values  of  a  are  chosen  to  follow 
percentages  of  the  peak  velocity,  as  shown  in  Figure  2.  Mathematically, 
the  transformation  is  as  follows: 
a'  *  a 


(6) 


8  m  f(z,t) 

where  a'  Is  the  depth  coordinate  of  the  new  coordinate  system. 

This  transformation  is  piecewise  linear  between  depths  where  velocity 
records  are  available.  Averaged  characteristic  information  available  from 
the  records  at  each  depth  is  used.  At  constant  depth,  the  time-like 
coordinate  a  is  related  to  t  through  a  series  of  straight  line  seqments  between 
a  and  t  values  established  for  particular  velocity  record  characteristics. 

Details  of  the  formulation  of  the  transformation  and  other  matters  concerning 
the  Interpolation  between  records  at  different  depths  are  given  in  the 
appendix. 

The  interpolation  of  the  velocity  between  depths  where  data  are 
available  Is  done  in  the  a',  a  coordinate  system  along  lines  of  constant 
a  values;  l.e.,  the  Interpolation  Is  performed  with  respect  to  specific 
velocity  record  characteristics.  Peak  velocities  at  adjacent  depths  where 
data  are  available  are  used  to  Interpolate  the  peak  velocity  at  In-between 
depths.  Thus,  half-peak  velocity  data  are  used  to  Interpolate  half-peak  velocities 
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and  so  on. 

In  order  to  establish  a  unique  definition  of  the  velocity  history  at  a 
given  depth,  the  velocity  record  Is  broken  up  into  a  series  of  segments  that 
begin  and  end  at  specific  a  coordinates.  The  segments  are  taken  to  be  the 
same  at  all  depths  forming  a  grid  work  in  the  z  \  a  space  as  shown  in 
Figure  3.  The  velocity  history  at  each  depth  and  over  each  segment  is 
represented  as  a  cubic  function  with  the  beginning  and  end  of  each  segment 
having  the  same  velocity  value  as  adjacent  segments  so  that  step  changes  in  velo¬ 
city  (infinite  acceleration)  are  ruled  out.  The  parameters  of  the  cubic  interpolation 
functions  are  evaluated  based  on  a  least  square  fit  to  the  velocity  records 
at  the  depth  in  question. 

Interpolation  between  depths  is  done  using  an  exponential  function 
that  begins  at  the  adjacent  upper  depth  and  ends  on  the  adjacent  lower  depth. 
Extrapolation  of  velocities  to  depths  above  the  shallowest  depth  for  which 
data  are  available  is  performed  by  extrapolating  the  exponential  interpolation 
function  derived  for  the  region  between  the  two  shallowest  depths.  Velocity 
at  depths  below  the  deepest  depth  for  which  data  are  available  is  obtained, 
in  a  similar  manner,  from  the  interpolation  function  of  the  two  deepest 
depths . 

Having  established  a  method  of  deriving  unique  velocity  records  for 
all  depths  and  times  from  the  measured  data,  stresses  and  strains  are 
calculated  from  Equations  (2)  and  (5)  using  central  differencing  and 
standard  numerical  Integration  techniques.  The  results  of  the  analysis 
have  been  found  to  be  rather  insensitive  to  the  discretization  used  in  the 
numerical  analysis.  The  major  factors  In  the  analysis  appear  to  be  the 
choice  of  velocity  histories  and  of  Interpolation  segments. 
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ANALYSIS  RESULTS 


Soil  stress  strain  relations  have  been  calculated  from  one-dlmenslonal 
finite  difference  calculation  results,  as  a  check  on  the  analysis,  and  from 
DISC  HEST  test  data.  The  finite  difference  calculation  considered  two  dry 
soil  layers  over  a  wet  soil  half-space  (Reference  1).  The  dry  soil  layers 
had  linear  loading  and  linear  unloading  moduli.  The  loading  modulus  for 
the  lower  dry  soil  material  was  1.5  times  as  stiff  as  for  the  surface 
layer;  the  unloading  moduli  in  both  layers  were  about  an  order  of  magnitude 
stiffer  than  the  loading  (Fig.  4).  The  velocity  histories  from  the  calcu¬ 
lation  were  Input  Into  the  analysis  developed  herein,  with  the  stress- 
strain  curve  results  shown  in  Figure  4.  For  this  analysis,  velocity 
histories  were  used  at  every  1.83  m  (6  ft)  of  depth  from  the  surface  down 
to  18.29  m  (60  ft)  and  every  3.05  m  (10  ft)  of  depth  from  18.29  m  (60  ft) 
down  to  45.72  m  (150  ft). 

The  technique  Is  able  to  track  properties  that  change  with  depth  and 
Is  able  to  follow  unloading  along  entirely  different  slopes  than  the 
loading  curve.  The  unloading  results  are  less  satisfactory  than  those  for 
the  loading  portion  of  the  response.  However,  because  the  unloading  is  so 
steep,  a  small  change  In  strain  can  make  a  large  difference  in  the  slope  of 
the  curve.  The  results  also  show  that  the  method  predicts  a  gradual,  rather 
than  a  sharp,  change  In  properties  with  depth.  This  Is  attributed  to 
smoothing  and  other  approximations  Inherent  In  the  use  of  exponential  and 
cubic  functions  to  perfbrm  the  required  Interpolations. 
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When  using  theoretical  results,  velocity  records  are  available  at 

* 

the  surface  as  well  as  at  large  depths.  With  experimental  data,  surface 
velocity  histories  are  not  known  and  histories  at  deep  depths  are  likely  to 
be  influenced  by  free  surface  reflections  from  outside  the  loaded  area 
that  produce  multi-dimensional  response  characteristics. 

After  gaining  experience  with  the  technique  using  analytical  results.  It 
was  then  applied  to  data  from  the  DISC  HEST  Test  I  event  (Reference  2) 
conducted  in  Ralston  Valley,  Nevada.  The  HEST  cavity  radius  on  this  test 
was  13.7  m  (45  ft),  and  data  was  obtained  at  eight  depths  down  to  15.2  m 
(50  ft)  (Fig.  5).  The  data  appeared  to  be  relatively  free  of  noise,  allowing 
most  of  the  records  to  be  Incorporated  in  the  evaluation  of  the  velocity  field 
history.  A  total  of  fourteen  records  over  the  eight  depths  (Fig.  6)  were 
Included  in  the  analysis.  Figure  7  shows  a  family  of  Interpolated  velocity 
histories  obtained  from  the  analysis. 

The  alrblast  pressure  history  was  measured  at  several  points  on  the 
ground  surface  within  the  HEST  cavity.  A  best  fit  through  the  pressure 
records  (Reference  2)  is  shown  in  Figure  8.  This  pressure  history 
was  Initially  used  in  Equation  (2)  to  calculate  stress  histories  at  depths 
of  Interest,  but  unsatisfactory  results  were  obtained.  The  main  problem 
was  that,  at  the  onset  of  Incipient  motion  (s  *  0)  at  some  depths,  the 
stress  was  nonzero  because  the  two  terms  on  the  right  hand  side  of 
Equation  (2)  did  not  exactly  cancel.  This  occurred  because  the  averaged 
surface  overpressure  and  the  velocity  field  interpolated  from  the  data  were 
not  completely  consistent. 

An  alternate  approach  is  to  compute  a  surface  pressure  loading  from 
the  velocity  field.  Setting  the  stress  equal  to  zero  in  Equation  (2) 
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at  the  arrival  of  motion  at  any  depth,  dQ  yields 

rdO(t) 

0  *  -p(t)+  I  p  dz 

fd0(*>  (7) 

p(t )  =  J  p  bz  dz 

A  pressure  history  derived  from  the  ground  motion  response  is  shown 
in  Figure  8  compared  with  the  best  fit  pressure  history.  The  pressure 
histories  agree  reasonably  well  in  the  early  time  of  the  motion.  The 
initial  slope  of  the  calculated  pressure  loading  is  not  as  steep  as  that 
obtained  from  measurements.  This  is  probably  due  to  smoothing  of  the 
velocity  histories  by  the  interpolation  process  and  the  difficulty  of  exactly 
predicting  ground  motion  at  and  near  the  surface  from  measurements  made 
below  the  surface.  The  impulse  histories  of  the  measured  and  calculated 
surface  loadings  are  also  shown  in  Figure  8.  They  compare  very  well  in  the 
early  time  of  the  motion  indicating  that  the  interpolation  process  averages 
out  variations  in  the  velocity  histories  in  a  manner  that  preserves  the 
overall  character  and  energy  content  of  the  response.  After  about  25  ms 
of  response,  the  calculated  and  measured  pressure  and  impulse  histories 
begin  a  significant  divergence.  This  can  be  attributed  to  the  reduction 
of  the  vertical  ground  accelerations  by  edge  effects.  The  divergence  of 
the  pressure  histories  may  serve  as  a  time  marker  for  the  demarcation  of 
when  two-dimensional  response  becomes  Important. 


Examination  of  the  data  in  Figure  7  shows  that  peak  velocity  is 
achieved  within  25  ms  down  to  a  depth  of  about  8  m.  For  larger  depths, 
the  time  to  peak  velocity  increases  with  depth  at  a  faster  rate  than 
might  be  expected  for  one-dimensional  response.  Information  on  properties 
at  deeper  depths  can  be  obtained  from  two-dimensional  finite  difference 
calculations,  however,  assumptions  are  required  to  obtain  the  uniaxial 
strain  properties  at  the  depths  because  the  response  is  two-dimensional  and 
essentially  only  vertical  motion  records  are  available. 

Figure  9  shows  the  interpolated  uniaxial  strain  stress-strain  plots 
obtained  from  the  velocity  field  and  the  derived  surface  overpressure. 

The  initial  slope  of  the  stress-strain  curves  and  the  4  MPa  secant  modulus 
are  plotted  as  a  function  of  depth  in  Figure  10  and  compared  with  the 
seismic  velocity  profile  of  the  site.  The  loading  properties  show  relatively 
small  variation  with  depth  down  to  9  m.  The  unloading  properties  are  not 
well  behaved,  but  most  of  the  unloading  occurs  after  the  25  ms  of  one¬ 
dimensional  simulation  time. 

The  initial  slope  of  the  interpolated  stress-strain  curves  produces  a 
modulus  that  compares  better  with  the  seismic  results  (except  very  near  the 
surface)  than  does  the  4  MPa  modulus.  This  is  to  be  expected  since  the 
moduli  obtained  from  seismic  measurements  are  representative  of  the  soil 
response  at  very  low  stress  values.  The  4  MPa  modulus  is  consistently 
lower  than  the  seismic  or  initial  slope  values.  This  reflects  the  softening 
of  the  soil  with  Increasing  stress,  characteristic  of  cemented  granular 
soil. 

The  seismic  profile  shows  a  soft  soil  layer  in  the  top  1.5  m  (5  ft) 
that  Is  not  present  In  the  results  of  the  interpolated  stress-strain  curves. 
It  Is  possible  that  the  material  Is  behaving  stlffer  than  would  be  expected 
from  the  seismic  profile,  because  of  strain  rate  effects.  However,  the 


9 


fact  that  the  Interpolated  peak  surface  pressure  is  lower  than  (about  15  percent) 
the  averaged  pressure  gage  data  indicates  that  the  near  surface  motions 
were  actually  larger  than  those  used  in  the  analysis,  which  would  result 
in  softer  near  surface  properties. 

Better  results  might  be  obtained  by  making  use  of  seismic  velocity 
information  in  extrapolating  velocity  field  data  to  obtain  surface  values. 

The  second  depth  that  velocity  data  was  recorded  in  DISC  HEST  Test  I  was  at 
1.5  m,  the  same  depth  that  a  sudden  hardening  of  soil  modulus  was  indicated 
by  seismic  data.  Thus  the  second  velocity  history  occured  at  a  transition 
region  in  soil  properties.  The  velocity  response  at  this  depth  then 
contained  information  that  was  more  characteristic  of  the  soil  below  the 
transition  boundary.  Therefore,  the  extrapolation  of  motion  field  to 
the  surface  was  influenced  by  the  second  seismic  layer.  An  alternate 
approach  to  performing  the  extrapolation  would  be  to  increase  the  peak 
surface  velocity  until  the  interpolated  peak  surface  pressure  was  in 
agreement  with  the  data.  In  cases  where  it  is  important  to  more  accurately 
define  the  material  properties  in  this  very  near  surface  layer,  then  gages 
at  two  depths  within  the  layer  (such  as  0.5  and  1.0  m)  should  be  used  in 
the  experiment. 
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CLOSURE 


The  purpose  of  this  paper  is  to  demonstrate  a  methodology  for 
determining  uniaxial  strain  mechanical  properties  of  soil  solely  from 
velocity  time  histories  obtained  from  a  high  explosive  field  test  event. 

It  is  similar  to  the  LASS  (Lagrangian  Analysis  of  Stress  and  Strain) 
methodology  developed  by  SRII  (Reference  3)  for  analysis  of  spherical 
motions.  However,  for  the  spherical  case  both  stress  and  velocity  data 
are  needed. 

Analysis  of  in  situ  field  test  data  is  generally  the  most  accurate 
technique  for  determining  In  situ  mechanical  properties.  The  material 
property  inversion  technique  described  herein  represents  a  first  step 
in  the  analysis  of  the  data;  the  complete  development  of  properties  at  a 
site  would  consider  all  available  relevant  Information,  such  as  seismic 
and  laboratory  test  results.  The  properties  derived  from  in  situ  data 
might  then  be  smoothed  or  adjusted  based  on  auxiliary  data,  as  long  as  these 
changes  were  within  the  uncertainties  of  the  in  situ  analysis.  These 
revised  properties  would  then  be  used  In  one-  and  two-dimensional  finite 
difference  calculations  to  verify  their  adequacy. 
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CONCLUSIONS 


(1)  Dynamic  stress-strain  properties  of  in  situ  soil  may  be  derived 
directly  from  velocity  histories  taken  from  surface  pressure  loading 

tests,  using  the  method  described  herein,  down  to  depths  where  the  ground  motion 
is  sufficiently  one-dimensional.  Results  have  been  obtained  from  experi¬ 
mental  data  for  the  DISC  HEST  Test  I  event. 

(2)  Surface  pressure  histories  may  be  derived  from  the  velocity  data  for 
the  duration  of  one-dimensional  response.  This  can  provide  a  check  on  the 
consistency  between  pressure  and  velocity  data.  The  time  when  the  measured 
and  derived  surface  pressure  loading  diverges  is  an  indication  of  the 
duration  of  one-dimensional  response.  In  the  DISC  HEST  Test  I  event,  the 
measured  and  derived  Impulse  histories  were  in  reasonable  agreement  out 

to  about  25  ms. 

(3)  The  Interpolation  functions  used  in'the  technique  are  effective  In 
deriving  soil  velocity  as  a  continuous  function  of  depth  and  time.  Since 
the  functions  are  based  on  measured  rather  than  hypothesized  soil  response 
characteristics  they  should  be  able  to  be  applied  to  tests  with  different 
soil  types  with  equal  success.  This  approach  can  also  be  generalized  to 
apply  to  motion  fields  that  are  dependent  on  two  spacial  coordinates. 
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Appendix  I  -  Derivation  of  Coordinate  Transformation  and  Interpolation 
Functions 

There  are  a  number  of  characteristic  times  that  are  clearly  Important  In 
describing  a  velocity  record  (Fig.  2).  The  two  times  of  dominant  Importance  are 
the  time  when  the  motion  begins  and  the  time  when  the  velocity  reaches  Its  greatest 
absolute  amplitude.  Other  points  such  as  the  end  of  the  record  and  time 
where  the  velocity  attains  given  fractions  of  the  peak  velocity  serve  to 
complete  a  listing  of  the  Important  characteristic  time  points  of  the 
record.  The  points  that  serve  to  best  describe  the  records  may  be  assigned 
labels  that  we  will  denote  by  the  symbol  s.  For  convenience,  the  labels 
can  be  made  numeric  and  assigned  values  that  Increase  with  time  for  a  given 
record.  By  convention  the  start  of  motion  is  at  8=0;  the  peak  velocity 
is  at  8=1.0;  and  the  end  of  the  record  is  assigned  a  large  number  such 
as  8*10.  Points  between  the  peak  velocity  and  the  end  of  the  record 
have  8  values  between  1  and  10  and  points  between  the  start  of  motion 
and  the  peak  velocity  have  8  values  between  0  and  1. 

Once  a  set  of  8  labels  have  been  chosen  they  can  be  applied  to  all 
of  the  records  at  all  of  the  depths.  The  value  of  a  remains  constant  with 
depth  along  a  characteristic  line  connecting  similar  time  points  of  different 
records.  At  a  given  depth,  points  with  particular  values  of  8  form  clusters. 

The  best  estimate  of  where  particular  a  points  ought  to  fall  can  be  obtained 
by  calculating  mean  values. 
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t(z.,8 )  »  \  ^  t  .(Z.,3 )  (A-l) 

v  ni  J  t 

•  • 

where  t(z_.,z)  Is  the  time  of  characteristic  point  a  at  depth  z.  and  record  j 
t(z.j8)  Is  the  average  time  of  characteristic  point  a  fbr  depth  a.;  and 
n.  Is  the  number  of  records  at  depth  z.. 

To  estimate  values  of  a  between  depths  where  records  exist,  a  straight 
line  can  be  drawn  between  the  average  time  of  the  various  a  labels  for 
existing  records.  In  this  way,  a  family  of  segmented  constant  a  curves 
may  be  obtained,  which  can  be  considered  as  a  new  time-depth  coordinate 
system.  Specific  features  of  the  velocity  response  history,  for  each 
depth,  occur  at  constant  values  of  the  time  like  a  coordinate.  The  symbols 
of  the  new  coordinates  af,a  are  related  to  the  ztt  coordinates  by  the 
following  transformation  relations. 

*'*3  (A-2) 

In  order  to  further  define  the  function  f(z,t)  the  variation  of  a  as 
a  function  of  t  at  a  given  depth  must  be  specified.  The  simplest  choice, 
and  the  one  that  will  be  shown  here,  Is  to  let  a  vary  as  a  linear  function 
of  t  between  each  of  the  labeled  values  of  a.  The  linear  function  has  the 
advantage  that  It  guarantees  that  a  unique  mapping  of  a  onto  t  coordinates 
exists  and  vice  versa.  This  linear  relation  between  a  and  t  Is  Illustrated 


IS 


In  Figure  A-1  and  shown  In  equation  form  below 
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for  l(Zla(i})  <t<  i(zt8(i+Vj 


where  a^  Is  the  value  of  s  at  the  t'th  characteristic  point  on  the 
velocity  history.  Equation  (A-3)  can  be  rearranged  to  give  the  solution 
of  t  as  a  function  of  e. 
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for  a(i)  <  a  <  a(i+1) 
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The  terms  t(z'v,a)  are  linear  functions  of  depth.  In  equation  form 


the  relation  for  t(zta^ )  Is 
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Equations  (A-3)  and  (A-5)  define  the  function  f(z,t)  of  equation  (A-2). 
.  The  Inverse  transformation  equation  Is 


t  -  g(z «; 

The  advantage  of  the  coordinate  system  Is  that  It  provides  a  convenient 
framework  for  Interpolating  velocities  between  depths.  Points  at  s-1 
will  be  related  to  peak  velocity  values  only;  points  with  a-1/2  will  be 
related  to  velocity  values  that  are  at  1/2  of  the  peak  value,  and  so  forth. 

At  each  depth  a  variety  of  records  are  typically  available.  For 
those  records  at  a  given  depth  that  are  valid,  we  have  the  problem  of 
forming  a  function  that  Is  representative  of  the  velocity  response 
history.  Since  the  records  are  complex  It  Is  Impractical  to  consider 
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using  one  equation  to  represent  the  entire  time  history  at  a  given 
depth.  An  equation  of  this  kind  would  likely  be  complicated  and  might 
vary  In  form  depending  on  the  depth.  A  more  realistic  approach  Is  to 
segment  the  velocity  record  and  treat  each  segment  Independently  of  the 
other  but  at  the  end  points  of  a  segment  have  velocity  values  that  match 
up  with  those  of  adjacent  segments.  For  convenience  we  will  require  that 
common  a  values  be  used  to  define  segment  boundaries.  It  Is  usually  better 
to  have  a  larger  number  of  segments  to  define  the  history  at  the  beginning 
of  the  motion  than  at  the  end  of  the  velocity  record.  The  time  history 
at  the  beginning  of  the  motion  changes  more  rapidly  than  at  the  end  and 
thus  should  be  more  carefully  described. 

Lagranglan  Interpolation  functions  are  particularly  well  suited  for 
establishing  the  velocity  histories  for  the  various  segments  of  the  response 
at  a  particular  depth.  These  kinds  of  functions  can  be  readily  defined  to 
any  desired  order,  but  the  higher  order  functions  have  unfavorable  properties. 
Since  the  function  must  describe  the  velocity  history  over  only  a  segment  of 
the  total  time  of  the  record,  a  cubic  function  should  be  adequate  to  give  a 
good  description  of  the  required  motion.  The  general  form  of  a  cubic 
Lagranglan  Interpolation  function  Is  as  follows: 

W  ,ai2-aU>  raU~ai2J 


(ail~ai3}  (ai2~ai3}  (ai4~ai2}  (ai (ai2~a i4}  (ai3~ai4) 


(A-7) 


Here  v^3  and  are  the  velocities  at  each  of  four  « 

coordinate  locations  a^t  a^,  and  a ^  within  segment  i  and  at  depth  j,  and 
Wfj  Is  the  function  defining  the  velocity  history  for  the  i'th  segment  and 
j'th  depth.  The  locations  a ^  and  a ^  will  be  considered  to  occur  at  the 


beginning  and  end  of  segment  i. 


At  the  point  where  s  Is  zero,  the  value  of  velocity  is  zero  also. 

This  is  always  true  by  definition  since  the  s«0  point  is  taken  to  be 

the  point  In  time  where  the  motion  response  begins.  Thus,  for  the  first 

Interpolation  segment  between  s*0  and  where  s2  Is  the  value  at  the 

end  of  the  first  segment,  the  value  v^1  Is  zero.  The  values  w„3 

and  v . ..  are  still  unknown. 
tj4 

For  any  given  time  history  at  a  particular  depth,  errors  will 
occur  between  the  velocity  given  by  Equation  (A-7)  and  the  time  history 
value.  For  a  given  value  a  the  error  can  be  written  as 


or 


*jkM  m  W8)  ‘  Wi/8) 


(A-8) 


E.y(8)  -  V.k(s)  -  hi2 -  \2ViQ2  ■  hisvijs  ~  hi4Vij4  (A'9) 


where  V^(a)  and  E^(a)  are  the  velocity  and  velocity  error,  respectively, 
for  time  history  k  at  coordinate  a  and  depth  j  and 
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The  error  value  S^(a)  may  be  either  positive  or  negative;  however, 
only  the  absolute  amplitude  of  the  error  value  Is  of  Importance.  The 
square  of  the  error  value  S^fa)  Is  always  positive  and  Is  then  a  better 


error  measure  for  the  purposes  of  judging  the  ability  of  the  function  w..(8) 

.  J'4 


to  fit  the  velocity-time  response  V..fa). 

v  * 

The  total  error  is  not  just  that  measured  at  one  coordinate  value 
for  one  record  but  is  the  integral  for  all  points  e  over  the  segment 


summed  for  all  of  the  records  available  at  a  given  depth, 
n . 
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We  must  select  velocities  v vijy  and  that  minimize  the  total 


error  E ^ We  will  assume  that  the  value  y^  is  always  constrained  by 
continuity  requirements  with  segment  i-1.  If  we  evaluate  the  segments  in 
order  starting  with  i*2,  we  need  only  match  the  constant  v^4  for  segment 
1  with  the  constant  y^  for  segment  2  to  preserve  continuity  of  velocity 
between  segment  1  and  2.  The  constant  is  always  zero  since  it  corresponds 
to  the  beginning  of  motion.  In  evaluating  the  constants  for  segment  1,  we 
need  only  find  values  for  Vjjz*  vlj3  and  vlj4'  ^  vij4  set  ec>ua*  t0 
then  for  the  second  segment  only  the  values  v2j2*  v2j3  and  v2j4  need 
be  evaluated.  The  process  can  be  continued  for  all  of  the  segments. 

The  positions  and  then  must  fall  at  the  beginning  and 
end  of  the  segment.  Thus  and  .  The  points  a^2  and 

can  be  chosen  as  equally  spaced  along  the  segment.  Thus 

9i-l 
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The  error  function  E ^  is  always  positive  and  Is  a  quadratic 
function  of  the  terms  vijs  an<*  vi*4  *  **  follows  that  its  minimum 

occurs  where  the  derivatives  of  E^  with  respect  to  u^2,  v^.-j  and  v^4 
are  all  zero.  These  conditions  produce  the  equations  for  the  determination 
of  the  unknown  velocity  constants. 
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*  hi3hi4Vij3  +  hi4Vij4^ds 
Evaluating  each  term  of  the  summation  and  integration  processes  of 

Equations  (A-14),  ( A- 15)  and  (A-16)  separately  and  placing  terms  of  known 

value  on  the  right  hand  side  of  the  equation  with  terms  of  unknown  value 

on  the  left  hand  side  of  the  equation  produces  the  following  results: 
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where 
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Equations  (A-17)  are  a  set  of  three  linear  relations  with  constant 
coefficients  and  three  unknowns  that  can  be  readily  solved.  The  solution 
is  not  a  major  problem  once  the  constant  coefficients  <?22,  <?23,  etc.,  and 
the  right  hand  side  constants  F^t  Fs  and  F  are  known.  However,  these 
constants  require  an  integration  that  is  not  trivial.  The  functions 
hi2*  hi3  and  hi4  are  not  slmPle»  and  the  velocity  functions  V are  typically 
defined  at  a  finite  number  of  points  rather  than  in  a  continuous  manner 
due  to  the  digital  definition  of  the  record.  The  integrations  for  the 
coefficients  of  the  equation  can  be  carried  out  numerically  with  no  difficulty 
since  the  functions  are  defined  for  all  values  a  within  the  segment.  A 
standard  procedure  for  performing  the  numerical  integration  is  to  divide 
up  the  segment  into  a  large  number  of  Intervals,  and  then  evaluate  the 
Integral  based  on  the  function  values  at  the  beginning  and  end  of  each  of 
the  Intervals. 

For  the  Integration  of  the  constants  F3  and  It  will  be  necessary 
to  define  values  of  rik  at  points  a  where  digitized  velocity  data  are  undefined. 
This  problem  Is  resolved  by  assuming  that  the  velocity  in  each  record  varies 


linearly  between  the  defined  values.  Figure  A-l  Illustrates  the  assumed 
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variation  of  the  velocity  record  with  respect  to  a  and  time.  By  making  the 
assumption  of  linearity  between  the  defined  velocity  record  points,  the 
velocity  record  becomes  in  effect  defined  at  all  points  and  the  integrations 
of  Equations  (A-24),  (A-25)  and  (A-26)  can  be  carried  out  to  obtain  the 
constants  F F^  and  F^. 

With  the  determination  of  the  constants  v ^  for  all  segnents  at 
all  depths,  a  set  of  velocity  records  are  available  at  each  depth  that 
are  representative  of  an  optimal  average  of  all  valid  velocity  records. 

The  problems  now  remaining  are  how  to  interpolate  between  depths  where 
velocity  records  are  recorded  and  also  how  to  interpolate  from  the 
velocity  record  at  the  shallowest  depth  up  to  the  surface  and  to  depths 
below  the  deepest  where  data  is  available. 

There  are  any  number  of  schemes  that  could  be  applied  to  the  inter¬ 
polation  of  the  velocity  records  between  depths.  C?;e  could  use  linear, 
quadratic,  cubic  or  higher  order  polynomials,  or  one  could  use  functions 
that  are  appealing  on  a  physical  basis.  Since  some  characteristics  of 
velocity  histories  are  known  to  decay  in  an  approximately  exponential  manner 
with  depth  within  a  given  material,  exponential  functions  should  provide 
useful  vehicles  for  interpolating  the  velocity  records.  After  some  experi¬ 
mentation,  it  was  found  that  the  exponential  functions  in  fact  produced 
more  favorable  results  than  the  polynomial  functions.  The  major  draw¬ 
back  to  the  polynomial  functions  is  their  tendency  to  oscillate.  This 
oscillation  produces  velocity  responses  that  can  Increase  with  respect 
to  depth  Instead  of  decrease  even  though  all  of  the  points  used  in  the 
Interpolation  show  a  decrease  in  velocity  with  respect  to  depth. 

The  simplest  type  of  exponential  fit  involves  placing  an  exponential 
function  between  two  points.  Suppose  it  is  required  to  determine  the 


variation  of  velocity  between  depths  j  and  along  a  time-like  coordinate 
line  s.  Using  Equation  (A-7),  the  velocities  w...  and  w..  are  obtained. 
The  equation  must  then  satisfy  the  following  constraints. 


(a+bz '.) 
-  e  J 


Wi(2+1) 


(A-27) 


The  constants  a  and  b  are  unknown  and  s',  and  s'.  ,  are  the  depth 

2  2+1  K 

values  at  depths  j  and  The  equation  can  be  solved  for  a  and  b  by 

taking  the  log  of  both  sides  of  the  equation. 


log  u . .  =  a  +  bz'. 

*  * 2  2 


log  v . ,  .  , ,  =  a  +  bz'.  . 
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Solving  for  a  and  b  produces  the  following 
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Knowing  the  constants  a  and  b,  the  velocity  w(z')  is  solved  from 
the  relation 


w(z')  -  e(a+hz’>  (z'...<z'<z\)  (A-31 ) 

2+ a  2 

A  special  problem  arises  when  the  velocities  are  negative  or  change 

sign  between  depths.  In  the  case  of  negative  velocities  w..  and  u...  . ,  in 

i(2+l) 

Equation  (A-29)  and  (A-30)  are  replaced  by  absolute  values  and  the  sign  of 
the  right  hand  side  of  Equation  (A-31)  Is  made  negative  Instead  of  positive. 
In  the  case  where  the  velocity  changes  sign  a  new  velocity  term  w'(z')  is 
defined  by  adding  or  subtracting  from  w(z')  a  constant  equal  to  twice  the 
absolute  value  of  the  velocity  at  the  lower  depth  point. 
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(A-32) 


w’Cz’)  -  w(z’)  +  2V}i(^+i)\ 

For  instance  1  f  w(z’ .)  is  negative  and  u(z'.  )  is  positive  then  2  jw.,.  7,(is 
subtracted  from  w(z')  to  compute  u'(z').  w(z’  .^)  would  then  be  equal  to 
-  wi ^  .  The  experimental  interpolation  is  computed  in  terms  of  w'(t’) 
which  always  has  the  same  sign,  w(z')  is  calculated  from  u'Ct'J  by  appro¬ 
priately  adding  or  subtracting  the  constant  from 


w(z')  »  u'fs';  ^  2  wi(j+i)  (A-33) 

The  selection  of  the  constant  is  arbitrary  but  is  not  a 

critical  one  since  sign  changes  in  velocity  occur  at  late  times  in  the 
velocity  histories  where  soil  motions  are  not  critical  in  determining  loading 
slopes  of  the  soil  stress-strain  behavior. 

To  extrapolate  velocities  to  the  surface,  the  velocity  records  of 
the  two  shallowest  depths  must  be  used.  The  only  piece  of  information 
available  at  the  surface  is  the  time  of  the  start  of  the  pressure  loading. 

The  first  s  label  line  at  e=0  can  then  be  drawn  from  the  first  depth  to 
this  point  at  the  surface.  The  s  label  at  s=l  should  fall  no  more  than 
a  millisecond  behind  the  label  for  8-0  since  the  rise  time  of  the  blast 
loading  is  very  short.  The  8  labels  falling  after  s»l  can  be  extended 
up  to  the  surface  based  on  their  slope  between  the  two  shallowest  depths. 
These  considerations  are  illustrated  in  Figure  (A-2). 

Having  established  the  8  coordinates  for  extrapolatinq  the  velocity. 

It  remains  to  establish  how  the  velocity  varies  in  the  upper  layer 
of  soil .  The  recommended  approach  is  to  simply  extend  the  exponential 
velocity  fit  between  the  two  shallowest  depths.  If  the  exponential  form 
of  the  function  has  a  physical  basis  then  this  form  ought  to  be  the  best 
possible  to  use  given  no  other  Information  about  the  velocity  In  this 
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region.  In  a  similar  manner  the  two  lowest  depths  with  data  are  used  to 
extrapolate  down  below  the  region  where  information  Is  available. 

The  velocity  time  history  for  all  times  and  depths  is  thus  obtained 
by  three  fundamental  steps.  First,  a  coordinate  transformation  between 
Zjt  and  z',8  coordinates  Is  established  where  certain  e  values  relate 
to  particular  features  of  the  velocity-time  record.  Second,  a  set  of  cubic 
Lagrangian  velocity  functions  are  established  at  each  depth  to  fit  a  set 
of  velocity  records  using  a  least  square  error  criteria.  Third,  an 
exponential  function  is  used  to  extrapolate  velocity  values  from  the 
Lagrangian  function  along  lines  of  constant  8  between  depths  where  the  velocity 
functions  are  defined. 
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Appendix  II  -  Notation 


The  following  symbols  are  used  in  this  paper: 


ail* ai2* ai3* ai4 


time-like  s  coordinate  locations  within  Interpolation 
segment  i. 


a,b  »  Unknown  constants  in  the  evaluation  of  the  exponential 
function  used  to  interpolate  between  depths. 


a..  ■  coefficients  of  set  of  linear  equations  used  in  least 
3  square  solution  for  y^.^. 


dQ(t)  *  shallowest  depth  at  time  t  where  soil  motion  and 
stress  due  to  soil  motion  are  zero. 


E.,  (s)  *  error  between  velocity  record  j  at  depth  k  and 
'  Lagrangian  interpolation  function  for  time-like 

coordinate  a. 


f(z,t)  *  function  defining  time  like  coordinate  s  in  terms 
of  a  and  t. 


g(z',a)  «  function  defining  time  t  in  terms  of  z'  and  a. 

*  comP°nents  °f  lagrangian  Interpolation  functions  u.Js). 

»  the  number  of  records  at  depth  i. 

p  •  pressure  at  surface 

e  ■  time-like  coordlrate  that  remains  constant  with  depth 
for  particular  features  of  the  velocity  response 
history. 

■  time-like  coordinates  at  the  beginning  and  end  of 
Interpolation  segment  i. 

t  •  time 

^(z.t8)  «  the  average  time  of  characteristic  point  a  for 
depth  a.. 

t.(z.,s)  •  the  time  of  characteristic  point  a  for  depth  z. 

3  and  record  j .  1 

u_  »  displacement  In  vertical  direction  (positive  direction 
Is  up). 

■  velocity  of  record  j  and  depth  k  for  time-like 

3  coordinate  e. 
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Appendix  I!  -  Notation  (continued) 


v.  •  velocity  in  vertical  direction  (positive  direction 
5  Is  down). 

v  • v . ■  velocity  values  used  with  Lagranglan  interpolation 
J  3  J  J  function  for  interpolation  segment  i  at  depth  j . 

»..(•)  •  Langrangian  interpolation  function  written  In  terms 
of  time-like  coordinate  0  for  interpolation  segment  i 
at  depth  j. 

z  *  vertical  coordinate  (positive  direction  is  down) 
coordinate  system]. 

0'  ■  vertical  coordinate  (z',9)  coordinate  system 
p  •  soil  density 

*  normal  stress  In  a  direction  (tensile  stress  Is  positive). 


27 


Appendix  III  -  References 


1.  Zelasko,  J.  S.,  "Sensitivity  of  MX-Relevant  Vertical  Ground  Shock 
Environments  to  Depth  to  Groundwater  and  Soil  Compressibility 
Variations",  U.S.  Army  Engineer  Waterways  Experiment  Station, 
April,  1980. 

2.  Jackson,  A.  E.,  et  al ,  "Ralston  Valley  Soil  Compressibility  Study 
Quick-Look  Report  for  DISC  Test  I",  U.S.  Army  Engineer  Waterways 
Experiment  Station  Structures  Laboratory,  Vicksburg,  Mississippi, 
May,  1981 . 

3.  Gradv,  D.  E.,  et  al,  "In  Situ  Constitutive  Relations  of  Soils 
and  Rocks",  DflA  3671Z,  SRI  International,  Menlo  Park,  California, 
March,  1974. 


DISCHEST  DATA 


-DISC  HEST  Test 


INTERPOLATION  /CHARACTERISTIC 

SEGMENT  \  /  LINE 


nvNiaaooo  .z 


,  s  Coordinate  Space  with  Interpolation  Segment  Boundaries 


VELOCITY  (m/s)  VELOCITY  (m/s)  VELOCITY  (m/s)  VELOCITY  lm/s) 


TIME  (ms) 


Fig.  6. -Velocity  Records  for  DISC  HEST  Test  I  Analysis 


TIME  (ms) 

Fig.  7. -Interpolated  Velocity  History  for  DISC  HEST  Test  I  Analysis 


WES  BEST  FIT  THROUGH 
SURFACE  PRESSURE 
DATA 


(DSS-BdW)  SSIDdWI 


(BdW)  3anSS3Bda3AO 


-Soil  Stress-Strain  Histories  for  DISC  HEST  Test 


0.001  sec 


V 


Label  Configuration  Between  Surface  and  Shallowest  Depths 


